Revisiting metric perturbations in tensor- vector-scalar theory 
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The cosmological behavior of modified gravity theories with additional degrees of freedom (DOFs) is typi- 
cally complex and can give rise to non-intuitive results. A possible way of exploring such theories is to consider 
appropriate parameterizations of these new DOFs. Here I suggest using the algebraic structure of trivial iden- 
tities, which typically occur at the level of the perturbed field equations, for defining such parameterizations. 
Choosing the example of Bekenstein's tensor- vector-scalar theory (TeVeS) and considering perturbations in the 
conformal Newtonian gauge, this parameterization is then used to study several aspects of the cosmological 
evolution in an Einstein-de Sitter universe. As a main result, I conclude that perturbations of the scalar field 
take a key role in generating enhanced growth if this enhancement is primarily associated with a gravitational 
slip. From this point of view, the previously found modified growth in TeVeS is truly a result of the complex 
interplay between both the scalar and the vector field. Since the occurrence of trivial identities of the above 
kind appears as a generic feature of modified gravity theories with extra DOFs, these parameterizations should 
generally prove useful to investigate the cosmological properties of other proposed modifications. As such pa- 
rameterizations capture the full nature of modifications by construction, they also provide a suitable framework 
for developing semi-analytic models of cosmologically interesting quantities like, for instance, the growth fac- 
tor, leading to various applications. Supplementary to numerical analysis, parameterizations based on trivial 
identities are thus an interesting tool to approach modified gravity theories with extra DOFs. 

PACS numbers: 04.50.Kd, 04.25.Nx, 98.80.-k 
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I. INTRODUCTION 

As is well known, the current concordance model of cos- 
mology has proven quite successful in reproducing observa- 
tions on the largest physical scales y]]. Based on general rel- 
ativity (GR), this model heavily relies on two phenomenolog- 
ically motivated ingredients, namely cold dark matter (CDM) 
and dark energy (DE), which are believed to account for ap- 
proximately 95% of the Universe's energy-matter content. 
However, there are still several problems. From a more fun- 
damental point of view, for instance, the introduction of DE 
remains theoretically challenging and appears extremely fine- 
tuned, despite the various proposals for its dynamics H. The 
CDM paradigm, on the other hand, struggles with several is- 
sues as well. Apart from the lack of direct experimental detec- 
tion yD, these include the question of its cosmological abun- 
dance 0] and problems related to the formation of structure 
on small scales J5|] which are still subject to intensive debates 
in the literature. 

Given this rather unsatisfactory situation, one might take 
a completely different view on the above. One possibility is 
to consider that the gravitational description genuinely differs 
from GR, mimicking the basic effects of the dark components 
and alleviating some, if not all, of the problems inherent to 
the standard framework. Additionally motivated by consid- 
erations beyond the field of cosmology, approaches into this 
direction have recently (re-)gained some interest in the scien- 
tific community and there exists now a plethora of modified 
gravity theories which could serve as an alternative to GR. 
Among many other proposals [6|, examples of such modifi- 
cations include conformal Weyl gravity Q, f(R) gravity Q8[| 
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or Brans-Dicke theory J9J]. In general, modifications of the 
gravitational sector come at the cost of simplicity, typically 
giving rise to complex and less transparent theories that ex- 
hibit a non-intuitive behavior and are difficult to work with. 
This is especially true if a modified framework introduces new 
independent degrees of freedom (DOFs), for instance, by con- 
sidering (non-minimal) couplings to additional fields such as 
scalars, vectors or higher-order tensors. Thus it is often hard 
to get a grasp on the detailed features and underlying mecha- 
nisms of such theories which are usually investigated with the 
help of numerical methods. Considering the limitations of a 
numerical treatment, i.e. the restriction to a finite set of theory 
parameters and initial conditions, it is only natural to look for 
alternative and supplementary ways of analyzing the impact 
of modifications on the cosmological evolution. A possible 
approach into this direction is to work out suitable parameter- 
izations of the new DOFs arising in a certain theory or gener- 
alized classes thereof. In particular, this might not only help 
to obtain a better understanding of a theory's intrinsic struc- 
ture, but also to highlight potential problems and limitations, 
which could provide a toehold for design improvements. It is 
almost needless to say that, for this purpose, such parameter- 
izations are required to be as close as possible to the actually 
assumed model. Note that this ansatz differs in motivation 
from the recent attempts of parameterizing deviations from 
GR in a model-independent fashion [10] which are interesting 
in their own right. 

Here I want to develop the above idea in the context 
of Bekenstein's tensor-vector-scalar theory (TeVeS) lUlll and 
use it to reinvestigate its behavior on cosmological scales. 
Considering perturbations around a spatially fiat Friedmann- 
Robertson-Walker (FRW) background, I will present an ap- 
proach toward parameterizing the new perturbation variables 
which is based on the occurrence of trivial identities at the 
level of the field equations. For the conformal Newtonian 



gauge, this parameterization is then used to study several 
aspects of the cosmological evolution such as the modified 
growth of density perturbations and the gravitational slip 
which have been identified as rather generic features of mod- 
ified gravity theories 111 211 . For simplicity and to allow an ana- 
lytic treatment of the background quantities, I will restrict the 
present analysis to an Einstein-de Sitter (EdS) universe which 
contains pressureless matter only. While this will suffice to 
gain a qualitative understanding of the perturbation evolution, 
it also provides an excellent approximation to the matter era 
of a realistic universe. 

In its first proposed form, TeVeS has been constructed 
as a relativistic extension for the modified Newtonian dy- 
namics (MOND) paradigm 11130 which aims at solving the 
missing mass problem on galactic scales by postulating an 
acceleration-dependent change of Newton's law. As such, 
TeVeS has been designed as a theory that is significantly dif- 
ferent from GR. Analyzing the rotation curves of spiral galax- 
ies, Milgrom estimated the characteristic scale of this transi- 
tion as flo ~ 1.2 x l(T 10 m s~ 2 , and it is interesting to note that 
MOND still appears to be in good agreement with much of 
the available extragalactic data [14]. As it turns out, the phe- 
nomenology associated with this kind of modifications and its 
generalizations H 1 511 is quite rich. It is therefore not too sur- 
prising that TeVeS and related theories have been subject to 
many different studies in the literature IU6I1 which range from 
the Solar System up to cosmological volumes. Note that, de- 
spite several inconsistencies arising from its construction, e.g. 
in the strong gravity regime II 1711 . the original formulation of 
TeVeS provides a viable theoretical framework for cosmolog- 
ical models ifflSiT - 

The paper is organized as follows: Starting with a brief re- 
view of TeVeS and its relation to MOND in Sec. HH the cor- 
responding cosmological background together with the modi- 
fied EdS model are discussed in Sec. [Hi] Considering pertur- 
bations in the conformal Newtonian gauge, I then introduce 
the parameterization of the new perturbation variables and in- 
vestigate the cosmological perturbation evolution in Sec. [TV] 
Finally, I conclude in Sec. [V] For clarity, some of the material 
involving lengthy expressions is presented in an appendix. To 
aid the reader through the technical parts, two tables of sym- 
bols are presented. Table ^summarizes all relevant quantities 
and parameters except for the ones related to perturbations 
which are separately listed in Table HH Throughout this work, 
I will assume a positive metric signature (+2) and use natural 
units, i.e. units where the speed of light, the reduced Planck 
constant and the Boltzmann constant are all equal to unity. If 
not stated otherwise, greek and latin indices run from to 3 
and from 1 to 3, respectively. 



II. FUNDAMENTALS OF TEVES 



Einstein metric g^ v , a time-like vector field A^ such that 

fV» = -L (!) 

and a scalar field <p. Furthermore, there is a second metric g^ 
which is needed for gravity-matter coupling only and obtained 
from the non-conformal relation 



guy = e 2 % v - lA^Av sinh(20). 



(2) 



The frames delineated by the metric fields g MV and g^ y will 
be called Einstein frame and matter frame, respectively. The 
geometric part of the action is exactly the same as in GR: 
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-J-g~d A x, 



(3) 



where R^, is the Ricci tensor of g MV and g the determinant of 
g MV . Note that the TeVeS constant G must not be mistaken for 
the Newtonian gravitational constant Gn (see Sec. Ill Bb . The 
vector field's action S v reads as follows: 



Sy 



with F„ 



1 



32ttG 



[[KsF^F^-AiA^ + l)] 



go x, 



(4) 



V^Ay - VyAfj and indices being raised and lowered 



with respect to g^, i.e. A M = g^Ay. Here the constant Kg 
describes the coupling of the vector field to gravity and A is a 
Lagrangian multiplier enforcing the normalization condition 
given by Eq. (Q]). Equation (01 corresponds to the classical 
Maxwell action, the field A^ now having an effective mass. 
The action S s of the scalar field involves an additional non- 
dynamical scalar field fi, and takes the form 



S s = - 



!6nG 



n 



niri^iy 



+ V(ju)] yRgd'x, 



(5) 



where W™ = g^ v - A^A V and VQi) is an initially arbitrary 
(potential) function. As the field /i is related to the invari- 
ant h^VjjipVyfi, however, it could in principle be eliminated 
from the action. Finally, matter is required to obey the weak 
equivalence principle, and thus the matter action is given by 



s m = f -C m [g,r B ,vr B ] ^d 4 x, 



(6) 



A. Fields and action 



where T B is a generic collection of matter fields. Note 
that world lines are by construction geodesies of the met- 
ric g^y rather than g^y. As usual, the corresponding equa- 
tions of motion can be derived by varying the total action 
S — S g + S „ + S s + S,„ with respect to the basic fields. De- 
spite its explicit bimetric construction, TeVeS may be written 
in pure tensor- vector form [20] and provides a particular ex- 
ample of generalized (extended) Aether-type theories I12l1l . 



Choice of the potential, quasistatic systems 
and relation to MOND 



TeVeS ifTll [1911 has originally been proposed as a bimetric 
theory of gravity which is based on three dynamical fields: an 



Whether or not TeVeS recovers the dynamics of MOND in 
the nonrelativistic limit depends on the assumed form of the 



potential V. Originally, Bekenstein made the choice 
3 ^o 



V(m) 



USnli 



[ft (4 + Ifi - 4£ 2 + jx 3 ) + 2 log (1 - A) 2 ] , (7) 



where the constant Ib corresponds to a length scale, p. - pi pa 
and yUo is a dimensionless constant 122(1 . One can show that the 
framework of GR is then recovered in the limit Kg — > and 
l B — > oo. To obtain the theory's nonrelativistic limit, one may 
apply the usual approximations for weak fields and quasistatic 
systems. In this case, one has V = dV/dp < 0, and therefore 
< p < fiQ. Using that also V{p) < for the given range, 
the resulting metric g„ v turns out to be basically identical to 
the metric obtained in GR if the nonrelativistic gravitational 
potential is replaced by 



77 _ „-2<pc 



(l+K B /2)- 



(8) 



where cpc is the cosmological value of <f> at the time the sys- 
tem in question breaks away from the cosmological expan- 
sion, and <D;v is the Newtonian potential generated by the mat- 
ter density p 12311 . In this approximation, it is consistent to 
assume that A^ is pointing into the direction of the timelike 
Killing vector associated with the static spacetime. Then we 
have 



hTVrfVyfi -» (V0) 2 = HV0H 2 . 
and the equation of the scalar field reduces to 
V • OuV0) = 8nGp. 



(9) 



(10) 



As has been shown in Ref. fi 111 , Eqs. ([H) and ( TTOb correspond 
to the MOND paradigm: If p — » po, the theory reaches its (ex- 
act) Newtonian limit, and the measured gravitational constant 
Gn is given by 



po + 2 - K B ( 
Ho(\-K B /2) 



(ID 



Similarly, the theory reaches its MONDian limit as p — > and 
the acceleration constant ao can be expressed in terms of the 
TeVeS and potential parameters, 



a = 



V6 e^ G 



(12) 



As can be seen from above, oo depends on 4>c an d mav there- 
fore, in principal, change with time. For viable cosmological 
models (see Sec. PH l. however, such changes are expected to 
be basically imperceptible 02411 . Note that requiring the emer- 
gence of MONDian dynamics fixes the potential only asymp- 
totically and leaves substantial freedom for constructing par- 
ticular models. 

It should be pointed out that potentials like the one speci- 
fied in Eq. (O exhibit a disconnection between the regimes 
relevant for cosmology and quasistatic systems, respectively. 
Since cosmological models require V' > fill] , one obtains 
// > jUq and thus cannot use the same potential branch as for 
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FIG. 1. Illustration of the generalized potential function V„(/j) given 
by Eq. J 1 3fc for n = 2 (solid line), 3 (dotted line), and 6 (dashed line). 



quasistatic systems (fi < /jo). Lacking a smooth transition be- 
tween these two regimes, however, it is unclear how bound 
systems such as galaxies would decouple from the Hubble 
flow or if such a decoupling results in the quasistatic limit 
discussed above. To resolve this issue, an interesting alterna- 
tive has been proposed in 02511 . with its cosmology studied in 
Ref. 12611 . In the present work, however, I will only focus on 
the theory's cosmological behavior and thus not follow this 
approach. 

Instead - for reasons that will become clear below - I will 
assume the following general class of potentials throughout 
this paper Il27ll : 



V n Qi) = 



Mi 



32jtI 



n + 4 + (n + \){1 



(n + l)(n + 2) 



(fi-2) 



n+\ 
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m=l 



(-D- 



(fi-2)" 



(13) 



where n > 2 IJ28I1 . Adopting different values of n, Fig. Q] 
illustrates the resulting potential shape as a function of fi. Note 
that the such generalized potential reduces to Bekenstein's toy 
model if n — 2. The derivative of V,,(/j) takes a simpler form 
and can be expressed as 



v'M 



3/4) „ 2 0" 

— fit- 
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1 



(14) 



Remembering that cosmological models must satisfy the con- 
dition V > and requiring that V is single-valued, one is al- 
ways free to choose between two possible pote ntial branches. 
In accordance with previous investigations 11111,1271 12911 . 1 shall 
use the branch ranging from the extremum at p — 2po to in- 



finity. Under these preliminaries, it was found that the poten- 
tial in Eq. ( [T3l gives rise to tracker solutions of the scalar 
field 1231, with a background evolution similar to other gen- 
eral cosmological theories involving tracker fields [30]. 1 will 
further elaborate on this behavior and an approximate analytic 
treatment in Sec. II11BI 



III. COSMOLOGICAL BACKGROUND 
A. Evolution equations 

Imposing the usual assumptions of an isotropic and homo- 
geneous spacetime, both g MV and g ^ are given by FRW met- 
rics with scale factors a and b — a exp <p, respectively, where <p 
is the background value of the scalar field ifllL Bill . Adopting 
a spatially flat universe, the modified Friedmann equation in 
the matter frame reads 



3// 2 = 87rG eff (p +p), 



(15) 



where the physical Hubble parameter is H — a/a 2 and the 
overdot denotes the derivative with respect to conformal time. 
Here p corresponds to the FRW background density of the 
fluid and the scalar field density takes the form 



(JJV' + V). 



r * 16nG 
The effective gravitational coupling strength is given by 

-2 



Teg 



Ge 



-40 



1 + 



d(f> 
dloga 



(16) 



(17) 



which is generally time-varying through its dependence on the 
scalar field </>. Just as in GR, the energy density p evolves 
according to 



p = — 3— (1 + w)p, 
a 



(18) 



where w is the equation-of-state (EoS) parameter of the fluid. 
In case of multiple background fluids, i.e. p = 2iPj, the rela- 
tive densities Q, are defined as 



Q, = 8nG 



eff; 



Pi 



Pi 



3H 2 p+p^ 
The evolution of the scalar field <p is governed by 



*=*&-*)- h 



377-0 + 4nGa 2 e- 4 * (p + 3P) 



(19) 



(20) 



where P is the fluid's background pressure, and the function 
U is related to the potential V, 

UQD^Jl + 2^. (21) 

In addition, the scalar field obeys the constraint equation 

f = I fl V 2 V (22) 



which can be inverted to obtain //(a, (/>, </>)■ For later use, I also 
introduce the relation 



a b 



a b 



_^2 

n4> 



4nGa 2 e' 4 



'(P + P) 



(23) 



which follows from combining Eq. (TT5l> wi th Eq. (1221 and the 
corresponding Raychaudhuri equation Bill. 

According to previous investigations lfllll27l 12911 . a broad 
range of expressions for the potential V, including the choice 
in Eq. (fT~3~T >. leads to exp <f> » 1 and p* <k 1 throughout cos- 
mological history. Therefore, the background evolution will 
be very similar to the standard case of GR, with only small 
corrections induced by the scalar field. 



B. Tracker solutions of the scalar field 

For the class of potentials specified in Eq. (fT~3~T >. it has been 
found that the scalar field exhibits a (stable) tracking behavior 
and synchronizes its energy density with the dominant com- 
ponent of the universe 1271 12911 . Tracking occurs as V tends 
to its zero point where /7 = 2|/o, and the evolution of the field 
<p during tracking is approximately given by 



= 00 + 



|l+3w| 



2p>o|l -w|-|l+3w| 



log a, 



(24) 



where (f> Q is an integration constant and/3 = ±1, with the ac- 
tual sign depending on the background fluid's EoS parameter 
w and Eq. d20b . Its density p# then exactly scales like that of 
the fluid, and the relative density parameter Q.^ turns approxi- 
mately into a constant, 



o* = 



(1 + 3wf 



6/4o (1 -w) 



2 - 



(25) 



Note that the right-hand side of Eq. ( 1241 slightly differs from 
the expression presented in Ref. 02711 . In Appendix [A] I dis- 
cuss why this is the case and show that Eq. (1241 is indeed the 
correct result. 

Following the lines of Ref. 12711 . Jl may then be expressed 
as p - 2yUo(l + e) with < e «c 1. Using V(2//o) = and 
expanding V to lowest order in e, Eq. d22l leads to 



e = 



1 (1674^2 



2 I 3/io a 



(26) 



It turns out that this is the only stage at which the constant lg 
enters the evolution equations. In preparation for Sec. [IV] I 
further take the time derivative of the above, which yields the 
useful relation 



-l 2/-2 -a -\ 
4>e=-U -4>-+cp)e. 
n \ a I 



(27) 



Note that stable tracking requires e to asymptotically decrease 
to zero, i.e. e — > 0. Therefore one has the condition e < 
which may be used to infer the proper sign of the parameter/? 
in Eq. (1241 (see Appendix lAli. 



IV. METRIC PERTURBATIONS IN TEVES 



A. Preliminaries 
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FIG. 2. Relative deviation of the Hubble expansion in the modified 
EdS cosmology from the ordinary GR case: Shown are the results 
foifio = 200 (dotted line), 500 (dashed line), and 1000 (solid line). 



C. Modified Einstein-de Sitter cosmology 

In what follows, I shall assume a universe entirely made of 
pressureless matter with perfect tracking of the scalar field, 
corresponding to the EdS model in GR. Setting P — w — 
fixes /? = —!, and thus the scalar field can be written as 



1 



2^o + l 



log a. 



(28) 



To find the proper value of ft, one may either insert Eq. 
into Eq. d20b . or use the argument presented in Appendix |A] 
Since the fluid evolves according to Eq. (fist , the density takes 
the form p = pga~ 3 , where p is the background density's 
value today. Thus exploiting Eq. ( 1251 ) allows one to rewrite 
the modified Friedmann equation in the matter frame as 



H 2 = // V 3+4/(2 "° 



where I have used the definition 



H~ = e~ 



, 87rGp 



1 + 



6/Uo 



1 



1 



2//o + l 



-2 



(29) 



(30) 



From Eq. d29) . it is evident that the deviation of the Hubble 
expansion from the ordinary EdS case is entirely characterized 
by the parameter pa. For several reasons pi ft. po should take 
a rather large value on the order of 100 - 1000, and thus this 
deviation will be very small. Assuming different choices of 
Po, Fig. |2]shows the relative difference between the models as 
a function of the scale factor a, indicating that the change of 
the expansion is only at the percent level. 



Now I will turn to metric perturbations around a spatially 
flat FRW spacetime in TeVeS. The starting point is the set 
of linear perturbation equations for TeVeS which have been 
derived in fully co variant form in Ref. Bill . For simplicity, I 
shall restrict the analysis to scalar modes only and work within 
the conformal Newtonian gauge. In this case, metric perturba- 
tions are characterized by two scalar potentials *F and <t>, and 
the line element in the matter frame is given by 

ds 2 = a 2 [-(1 + 2T)c/t 2 + (1 - 2<I>)£y<foW] . (31) 

Similarly, one needs to consider perturbations of the other 
fields: While the fluid perturbation variables are defined in 
the usual way, i.e. the density perturbation, for instance, is 
expressed in terms of the density contrast 6, 



p=p + 5p=p(l+5), 
the scalar field is perturbed as 

(p = cp + ip, 



(32) 



(33) 



where tp is the scalar field perturbation. Finally, the perturbed 
vector field is written as 



A,, = ae * [A^ + cty) , 



where A„ = (1,0,0,0) and 



a M = Q¥- ip, Va) . 



(34) 



(35) 



Note that the time component of the vector field perturbation 
is constrained to be a combination of metric and scalar field 
perturbations, which is a consequence of the unit-norm con- 
dition in Eq. ([T). Therefore, one needs to consider only the 
longitudinal perturbation component a. 



2. Einstein-frame perturbations 

Instead of using Eq. QTT l, one may also express pertur- 
bations in the Einstein frame 11291 [3111 . In this case, metric 
perturbations are written as 



|oo = -fcV 4 ^(l+2>P), 


(36) 


got = -b 2 dil 


(37) 


gy = ft 2 (l-2*)*y. 


(38) 



In terms of matter-frame variables, the Einstein-frame pertur- 
bations are given by 



Y = ¥ - tp, 

(J) = O - <p t 



(39) 
(40) 

(41) 



To avoid lengthy expressions in the perturbed field equations, 
it is convenient to work with variables from both frames. As I 
will frequently use these equations in the following sections, 
the full set of linear perturbation equations in the conformal 
Newtonian gauge is, for clarity, given in Appendix 151 

Since the equation governing the evolution of the scalar 
field perturbation is of second order, it is further helpful to 
introduce an auxiliary field y which allows one to split the 
scalar field equation into a system of two first-order equations 
Oil . Performing this split for the present gauge choice, the 
field y is given by the relation (see Appendix iBli 



y 



rfab*), 



(42) 



where the function T has the form specified in Eq. (!B4b and 
the g, denote a collection of background quantities and their 
derivatives with respect to conformal time. 



B. Parameterizing the new degrees of freedom 

To begin with, one may take a naive approach and try to 
relate the new perturbation variables to perturbations of the 
metric or the matter fluid by a suitable combination of the 
field equations. Considering the perturbation equations in Ap- 
pendix [B] and assuming a general matter fluid whose back- 
ground evolution is given by Eq. (fTST l, I take the time deriva- 
tive of Eq. ( IB8I 1 and eliminate 6 with the help of Eq. (IB2b 
from the resulting expression. The next step is to get rid of 
the time derivatives of <t>. This can be achieved by exploiting 
an algebraic relation which is obtained from combining Eqs. 
( |B91 l and JBlOt . Finally, using Eqs. d39b and dB4b . one arrives 
at the equation 



abb —2 - .-j 

2 Jt<p - 47rGaV 40 p(l + w) 

a b b 



¥ 



JiU(f> + 20- + A + fi(f> + 47rGaV 40 p(l + 3w) 



(43) 



<P- 



From Eq. d27] i, one immediately sees that the coefficient in 
front of "T vanishes. Thus the above gives a trivial identity 
and one cannot infer any information on the relation between 
the scalar field perturbation <p and the metric potential W. This 
somewhat reflects the fact that if corresponds to a full DOF 
in the theory and the occurrence of trivial relations like in Eq. 
d43l is indeed a generic feature of modified gravity theories 
which introduce additional DOFs Q2I1 . 

The idea is now to use the algebraic structure of these trivial 
identities as a starting point for defining a suitable parameter- 
ization. For purposes that will become clear below, consider a 



new function B l0 such that 



¥ = B v ip. 



(44) 



In general, B^ will change with time and exhibit scale- 
dependence which, for instance, may come as a ratio to some 
preferred scale(s) of the theory. Moreover, it will also depend 
on the used cosmological background, the theory parameters, 
the scalar potential V(fi), and the choice of initial conditions. 



As for the perturbation variables y and X, it is possible to 
motivate a relation similar to the above. To see this, I multiply 
Eq. (IB7b with the scale factor b and take the time derivative. 
Combining the result with Eqs. (IB5b and ( IB 10b , I eliminate E 
and the time derivative of (,, respectively. Substituting (p and y 
with help of Eqs. (IB3b and (1B4I) . respectively, one eventually 
ends up with 



b 11 — a— -J 
3 £0__0 + 

b U a 



■ cf> -\ 4nGa z e 



2„-4<4- 



p(l+3w) 



Ji4> + Ji(f> + 2 j 4> + - j JI<p + 4;rGa V 4 ^p( 1 + 3w) 



-2 



.b U 



- 1 + - + JM + 47rGaV 4 ^p(l + w) 
bj b 

(45) 

where k is the conformal wave vector and k = |k|, Using the 
background relations presented in Sec. IIII Al one finds that 
Eq. d45b again yields a trivial identity, with a form very sim- 
ilar to that found before. A direct comparison between Eqs. 
d43b and d45T > suggests the definition of another function B y 
which relates y to % and the time derivative of <$>. A suitable 
expression may be obtained by formally subtracting Eq. (l44l 
from Eq. (|43T > and "normalizing" the resulting equation, lead- 
ing to the ansatz 



3$ + k l l + B y y = 0. 



(46) 



The expressions given by Eqs. d44t and (l46l i provide formal 
closure relations which may be used to eliminate perturbations 
of the scalar and vector field from the evolution equations (see 
Appendix IClfor details). Thus perturbations to the new DOFs 
are genuinely absorbed into the functions B^ and B y . While 
the function B v simply describes the magnitude of scalar field 
perturbations relative to the metric potential W, the interpre- 
tation of By is less obvious. As the term k~l appears in its 
definition, however, one may expect that B y characterizes the 
impact of vector field perturbations. As will become clear in 
the following sections, this is to some extent the case and the 
current parameterization disentangles the effects of the vec- 
tor and scalar fields. Although the functions B^ and B y are 
entirely determined by choosing a specific model and a set of 
initial conditions, one may view them as "free parameters" de- 
scribing modifications induced by the additional DOFs. Con- 
cerning the analysis presented below, I will assume that B^ 
and B y are regular and require that B^, B y + at all cosmo- 
logical times and scales, which is necessary to exclude patho- 
logical cases. 

Although I have only considered the situation within the 
conformal Newtonian gauge, note that the occurrence of triv- 
ial identities is gauge-independent, and thus there is always 
the possibility of constructing similar parameterizations for 
different gauge choices. As parameterizations of this kind 
could prove useful for other modified theories, it might fur- 
ther be interesting to adopt a gauge-invariant approach toward 
such parameterizations, which I leave for future work. 



C. Subhorizon scales 

In the following, I shall assume the previously discussed 
modified EdS cosmology with perfect tracking of the scalar 
field. This allows one to use the corresponding background 
expressions presented in Sec. IHICI and considerably simpli- 
fies the analysis of the modified equations. Since /7 resides 
close to its minimum in this case, i.e. Ji = 2yUo(l + e) with 
e« 1, one may further exploit the two first-order expressions 
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l=l- 2 -e 
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1 e, 
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(48) 



which is useful to rewrite terms involving the field U. 



1 . Modified potentials 




Adopting an EdS universe together with the closure rela- 
tions presented in Sec. IIVBI one may now write metric per- 
turbations solely in terms of the matter fluid variables, B^ and 
By. The resulting expressions are quite lengthy and can be 
found in Appendix|C] As a first application, I shall investigate 
the theory's behavior on scales much smaller than the horizon 
where aH/k «; 1 . To allow further progress, I will addition- 
ally assume that the typical time variation of B^ and B y is 
comparable or smaller than that of the Hubble expansion, i.e. 



<aH, 



<aH, 



(49) 



such that the contribution of \B l/y /BJ and \B y /B y \ compared 
to k is much smaller than unity II 3311 . Inserting the logarith- 
mic approximation for <p specified by Eq. d28l l and using the 
assumptions given in Eq. d49l , I expand the corresponding 
equations for the Einstein-frame potentials T and €> in powers 
of aH/k. To lowest order, this yields 
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where 



A=- 



e 4 % [2 (B y + 2/iq) (e 4 * - l) + K B B y ] 
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FIG. 3. Evolution of subhorizon density perturbations in the modi- 
fied EdS cosmology: Assuming the model parameters specified in 
the text with B v = B r = 3/jo/2, the figure illustrates results for 
K B = 1 (dashed line), 0.1 (dotted line), and 0.09 (dashed-dotted 
line). For comparison, the corresponding evolution in the ordinary 
EdS model is also shown (solid line). 



Note that although I have omitted the resulting expressions to 
first order in e for clarity, their contribution, which typically 
amounts to a correction of a few percent or less, will fully be 
taken into account when conducting calculations. 

Compared to the situation in GR, where A = B — 3/2 for 
the EdS model, the above immediately shows that the poten- 
tials exhibit a more sophisticated time dependence. This de- 
pendence does not only involve the perturbation quantities B^ 
and B y , but also the background scalar field <p which is there- 
fore expected to have a significant impact on the potentials' 
evolution and the growth of density perturbations. 



2. Growth of density perturbations 

Equipped with an analytic expression for the the potential T* 
(or equivalently *?), I now proceed with the analysis of struc- 
ture growth in the context of TeVeS. As is well known, the 
ordinary EdS model in GR gives rise to a growth equation of 
the form 



(fid 3_dS_ _3_ 
da 2 2a da 2a 2 



(54) 



with the two solutions S oc a~ 3 ^ 2 and 6 oc a. Following the 
same derivation as in GR, the TeVeS analog of Eq. ( |54l for 
the present assumptions reads 

d 2 5 1 / 4 \ d8 A , A 



where A is given by Eq. d52l and <p evolves according to Eq. 
( l28l l. It is instructive to further simplify the expression for A 
by exploiting that the background field (p is much smaller than 
unity, i.e. \<p\ «: 1. A straightforward expansion immediately 
yields 
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(56) 



where I have additionally neglected terms proportional to e 
and used that fig » 1. As is obvious from the above, the 
vector field enters the growth equation not only through B y , 
but also implicitly through Kb- 

Adopting the potential specified in Eq. ©, i.e. the n — 2 
case, numerical studies of the Te VeS cosmology revealed that 
structure is allowed to form more efficiently than in the frame- 
work of GR 11290 . A later analysis of this behavior has identi- 
fied vector perturbations and their associated instability as the 
key ingredient for the found enhanced growth, where it has 
further been shown that this enhancement only occurs if the 
coupling constant Kb is chosen small enough 041 [3511 . On the 
other hand, it has been argued that scalar field perturbations 
play only a negligible role for structure formation, with ip os- 
cillating around a value roughly fio times smaller than metric 
perturbations during the matter era. Here I seek a detailed un- 
derstanding of this enhanced growth mechanism based on the 
modified growth equation and the functions B^ and B y . 

In accordance with the above findings, I will assume from 
here on that B^ takes values on the order of po, i.e. B^ ~ fio. 
To start with a simple case, consider first the situation that A 
is dominated by the zeroth-order term such that A » 3B^/2/jq, 
and that B^ « const. Then Eq. ( |55] l can be solved analytically 
and the growing mode approximately evolves as a power law, 
6 oc a p with 



1 
P = 4 



1 + 24^ - 
1 j"o 



\ 



(57) 



where the ordinary EdS behavior is obviously recovered if 
B^ - yU<). Next, one may ask how the first-order contribu- 
tions to A modify the growth of density perturbations. To an- 
swer this question, one first observes that, since B^ ~ jiq and 
|0| <K 1, the relevant terms in Eq. (|56T > are those depending 
on the function B y . Cleary, these terms can have a significant 
impact on A only if either B y or the product KbB 7 acquire suf- 
ficiently small values which are able to cancel the suppression 
by <f>. The resulting modification, i.e. an effective increase or 
decrease in growth, then depends on the actual numbers that 
are (or need to be) plugged in. 

To further simplify the present analysis, assume that any 
relevant first-order contribution in Eq. (1561 comes solely from 
the term proportional to K B l . This property may formally be 
implemented by imposing the additional condition \B y \ > B v . 
As will become clear in the following section, such an as- 
sumption is not unreasonable if the modification of growth is 
associated with a gravitational slip, i.e. a significant difference 
between the two metric potentials T* and <I>. For the same rea- 
son (see Sec. IIVC3I ), I will here exclude the case B^ - fio 
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FIG. 4. Same as Fig. [5] but now assuming a configuration with 



B 9 = 9yu /10 and B y 



-9^o/5: Shown are results for K B = 1 



(dashed line), 0.005 (dotted line), and 0.001 (dashed-dotted line), 
together with the ordinary EdS case (solid line). 



for which the first-order term under consideration vanishes. 
Then, for a fixed sign of and B^ + /iq, the form of Eq. (l56i l 
suggests that there are two principle configurations which can 
lead to an increase of growth. Choosing <p < 0, for instance, 
one finds that either B^ > [1q and B y > or B^ < /uo and 
B y < are required to obtain the desired enhancement (Note 
that a small negative value for the field cp is in harmony with 
the results of Ref. 0411 and, contrary to previous claims, does 
not automatically violate causality [36]). 

To demonstrate that this behavior truly emerges from the 
full potential expression in Eq. ( 1521 . one may numerically 
solve Eq. d55l l for different sets of trial parameters. Since the 
aim here is to develop a qualitative understanding, I will con- 
sider the case that both B^ and B y are approximately constant 
M37I1 . Assuming a B y > scenario with B ip -B y — 3jt/o/2 and 
setting ju = 1000, l B = lOOMpc, O = -0.003, the Hubble 
constant Ho = 100 km/s/Mpc, and n — 4 for the scalar poten- 
tial, Fig. [3] illustrates the numerically calculated evolution of 
5 for different values of Kb and an arbitrary, but fixed choice 
of initial conditions at a — 0.01. As can be seen from the 
figure, this simplistic model recovers an increased growth for 
small values of Kb (< 0.1). For large values of Kb (> 1), how- 
ever, this enhancement does not occur and the density con- 
trast follows a power law, with an index p given by Eq. 07] ). 
Similar results can be obtained for B y < 0, and leaving the 
other parameters unchanged, an example for B^ = 9/j.q/IO 
and B y - -IB^ - -9fio/5 is shown in Fig. [4] 

While both configurations, i.e. B y > or B y < 0, can 
give rise to increased growth, they substantially differ in the 
way this growth is sourced by the potential *¥, If B y > 0, 
models of the above kind suggest that the potential ¥ evolves 
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FIG. 5. Gravitational slip associated with the simplistic growth 
model shown in Fig. [4] The figure illustrates results for K B = 1 (solid 
line), 0.01 (dashed line), 0.005 (dotted line), and 0.001 (dashed- 
dotted line). 



in an exponential fashion. This, in turn, leads to a faster and 
faster increase of the density contrast as seen in Fig. [3] What 
is disturbing is that all such models with <p < eventually 
run into a singularity, which appears as a consequence of the 
denominator in Eq. ( |52| |. Although one may imagine that 
realistic time-varying choices of B^ and B y avoid a singularity, 
this result provides a strong indication against the B y > 
case. On the other hand, models with B y < are generally 
free of this problem and associated with a slower, logarithmic 
growth of the potential. Using the gravitational slip in the 
following section, it will turn out that the enhanced growth 
reported in Ref. [34]] does indeed correspond to the situation 
where B v < 0. 



3. Gravitational slip 

To conclude the discussion of subhorizon scales, I will 
briefly demonstrate how the previously discussed mechanism 
giving rise to enhanced growth generates differences between 
the matter-frame potentials M^ and O. Remember that in GR, 
such a difference can only be caused by anisotropic stress of 
the matter fluid which is basically negligible for the cosmo- 
logical evolution at late times. For this purpose, it is useful to 
introduce a quantity ^q characterizing the difference between 
the two metric potentials. Here I choose the definition 
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Assuming that B^ ~ jxq » 1 and 
the above expression becomes 



«: 1 (see Sec. IIVC2I ), 
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where I have again neglected terms of order e. Interestingly, 
the first-order contribution in Eq. (|59l does not depend on the 
function B y and contains a term proportional to K^iB^ - fio), 
strongly resembling that previously found in the expansion of 
the function A. Now consider the case B^ = jj.q for which 
this contribution disappears. From Eqs. (l52t and (l5Jt , it then 
follows that the zeroth-order expressions for A and B are al- 
most identical, only differing by a factor of exp (40) which 
appears in the second term of the numerator. Consequently, 
the relative difference between the metric potentials will only 
involve terms on the order of (p and e, becoming vanishingly 
small independent of what is assumed for B y . Therefore, if the 
enhanced growth is primarily sourced or driven by the gravi- 
tational slip, this heavily suggests that one must have B^ + jUo 
and picks out the term oc K B l in Eq. ( f5o*b as the predominant 
contribution, in accordance with the growth analysis presented 
in the previous section. 

Since the sign of ^q eventually depends on whether B^ takes 
values that are smaller or larger than fio, one may further use 
the gravitational slip to infer the sign of B y and thus the proper 
growth configuration from the numerical results obtained in 
Ref. 0411 . There it has been found that t;c significantly grows 
during the matter era with £c > and <p < Q8J1 . If one 
uses Eq. d59l , this leads to B^ < hq and thus B y < as 
increased growth would not occur otherwise (again, see Sec. 
IIVC 21 ). As such, the additional growth previously found in 
Ref. 112911 generally corresponds to the situation of the sim- 
plistic model depicted in Fig. @] Assuming different values of 
Kb, the slip parameter ^q associated with this particular model 
is illustrated in Fig. [5] 

From the findings of Sec. IIVC 21 and the above discussion, 
it becomes clear that the modified growth in TeVeS is truly a 
result of the complex interplay between the newly introduced 
DOFs, i.e. both the scalar and the vector field 0911 . While 
the scalar field enters the growth equation through B^ and its 
background value (f>, the vector field affects the evolution of 
6 through perturbations, i.e. the function B y , and the cou- 
pling constant Kb- If the observed enhanced growth is mostly 
sourced by the gravitational slip, perturbations of the scalar 
field take a key role in the growth mechanism and need to 
evolve in a way such that B v takes values sufficiently different 
from//o- 



D. Superhorizon scales 

Similar to the previous section, it is possible to analyze the 
theory on scales much larger than the horizon. In this case, 
terms proportional to k 2 may safely be neglected in the pertur- 
bation equations, i.e. k — > 0. Using Eqs. ( IB7b and ( IB lb . one 
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thus starts from 



„b U b-\ ae-*- 



47rGflV 40 p(5-2 V ) (60) 



and 



5 = 30 



3 (6 + ip) . 



(61) 



As usual, the above implies that the combination 6 - 3<5 is 
conserved over time. Using Eqs. (l60t and (|6TT > together with 
the remaining perturbation equations and the closure relations 
from Sec. IIVBI it is possible to arrive at an equation govern- 
ing the evolution of the Einstein-frame potential O. The corre- 
sponding derivation is sketched in Appendix IC 21 and the final 
expression takes the form 



CD + Co-<t> = 0, 
a 



(62) 



where the coefficient Co depends on background quantities, 
the parameter functions B v and B y , and their time derivatives 
(again, see Appendix IC 2\ . Since there appear only derivatives 
of the metric potential <t> in Eq. j62l , one immediately obtains 
the solution O = const. The second solution is determined by 
the evolution of Co in the modified EdS cosmology, but find- 
ing it requires detailed knowledge about the functional form 
of both B v and B y . To get some more insight, however, con- 
sider the situation that the time derivatives of these functions 
are approximately negligible. In this case, one can show that 
the coefficient Co approximately turns into 



Co 



By — IBy + 1 



2ZL 



B y + 2' 



(63) 



where I have again assumed that jUq » 1 and e «: 1 . The 
corresponding solution is then given by a power law, <J> oc a 1 ', 
with index 



5v 



2 B y -2{B ip + l) 



(64) 



The nature of this solution is obviously fixed by the actual 
configuration of the functions B^ and B y . To end up with a 
decaying mode, for instance, Eq. (l64l yields the condition 
2{B ip +l)/B y < 1. Remembering that scalar field perturbations 
should appear suppressed compared to the metric potentials, 
i.e. \B lf \ » 1, this is typically satisfied if B^ and B y are of 
opposite sign, but there is also substantial freedom for other 
configurations. In particular, Eq. (|64l reduces to the ordinary 
GR case if B y - 5(B^ + 1). Taking the view that TeVeS mod- 
ifies GR without introducing any very radical changes in the 
principal behavior of this solution, it seems reasonable to as- 
sume that it will generally correspond to a decaying mode for 
viable and realistic choices of Co- If this holds true, one is left 
with the result that the Einstein-frame potential <t> is frozen. 

As can be seen from Eq. (16 It . however, this result does 
not necessarily apply to the matter-frame potential <I> and the 
evolution of the density contrast which additionally involve 
the scalar field perturbation tp. Following a similar procedure 



TABLE I. Summary of relevant generic symbols and parameters: 
Note that barred symbols denote the corresponding background 
quantities in the FRW spacetime, and an additional subscript "0" in- 
dicates the value of a background quantity today (a =1). 



Symbol 



Description 



g m Metric tensor in the matter frame 

g m Metric tensor in the Einstein frame 

a Scale factor in the matter frame 

b Scale factor in the Einstein frame, b = ae* 

<p TeVeS scalar field 

Af TeVeS vector field 

F^, Field strength tensor of the TeVeS vector field A 11 

p Auxiliary nondynamical scalar field 

VQS) Free potential function of the scalar field p 

U(p) Function of p related to V, U = p + 2V/V" 

p Density of the matter fluid 

P Pressure of the matter fluid 

w Matter fluid's equation-of-state parameter, w = P/p 

H Physical Hubble parameter 

G Bare gravitational constant 

G e ff Effective gravitational coupling for FRW dynamics 

K B Coupling constant of the TeVeS vector field A^ 1 

V n (p) Generalized Bekenstein potential, see Eq. J 1 3b 

n Integer characterizing the potential V„,n > 2 

p Dimensionless parameter of the potential V„ 

I B Overall length scale entering the potential V„ 

e Offset of Jl/2p from minimum during tracking, e <k 1 



as for Eq. (I62t . it is also possible to obtain an equation for ip 
which takes the form of a simple mechanical oscillator, i.e. 



ip + d\(p + doip — 0. 



(65) 



Inserting the background relations of the modified EdS cos- 
mology and using fio » 1, the coefficients in Eq. d65l l can be 
expressed as 



2 + 



B, 

2/io/ a 



B,p a 

2pQ a 



(66) 



to lowest order in e. In general, one has d\ , di + 0, and thus 
solutions to Eq. j65l ) will vary over time. Therefore, mak- 
ing use of Eq. drjTT l, the potential <t> and 5 are not strictly 
conserved on superhorizon scales. The actual behavior of the 
perturbation ip is determined by the quantity d\ - 4-do, leading 
to damped oscillations or exponentially growing and decaying 
solutions. Given the typically found magnitudes of scalar field 
perturbations Il29ll34l 13511 . however, one may expect a rather 
small impact on the evolution of <$> and 6, but this is generally 
subject to initial conditions. Finally, note that such solutions 
with ip + also indicate an evolution of the Einstein-frame 
potential V F. This can be seen from 



d 1 d ~ 

ip + 2-<p+ T^O, 

a 2p.Q a 

which follows from combining Eqs. (1441) . (l65T l and (1551 . 



(67) 
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V. CONCLUSIONS 

The cosmological behavior of modified gravity theories 
with extra DOFs is typically complex and can give rise to non- 
intuitive results. A possible approach toward exploring such 
theories and their modifications is to consider suitable param- 
eterizations of these new DOFs. In this work, I have suggested 
using the algebraic structure of trivial identities, which typi- 
cally occur at the level of the perturbed field equations, as a 
starting point for defining such parameterizations. Choosing 
the specific example of Bekenstein's TeVeS theory and con- 
sidering only scalar modes in the conformal Newtonian gauge, 
I have used this ansatz to construct formal closure relations 
for the additional perturbation variables arising from the extra 
DOFs. These relations introduce two new auxiliary pertur- 
bation quantities, the parameter functions B^ and B y , which 
allow one to re-express perturbations of the TeVeS scalar and 
vector fields in the evolution equations and prove useful when 
investigating different cosmological aspects of the theory. 

Assuming the modified EdS background cosmology, i.e. 
the background solution of a universe dominated by pressure- 
less matter together with the cosmological scalar field acting 
as a perfect tracker field, I have studied in detail the theory's 
behavior on subhorizon scales. Focusing on the evolution of 
density perturbations, the present parameterization has been 
adopted to derive the TeVeS analog of the well-known growth 
equation and to address the question of how the additional 
DOFs affect the dynamics of the density contrast 6. Identi- 
fying the relevant terms of the TeVeS growth mechanism, I 
have explicitly demonstrated which configurations can lead 
to an increased growth rate and discussed their connection 
to the gravitational slip, i.e. a significant difference between 
the two metric potentials. As a main result, I conclude that 
perturbations of the scalar field take a key role in generating 
enhanced growth if this enhancement is primarily associated 
with a gravitational slip. From this point of view, the previ- 
ously found modified growth in TeVeS is truly a result of the 
complex interplay between both the scalar and the vector field. 
Considering the limit of superhorizon scales, it can be argued 
that non-decaying solutions of the Einstein-frame potential O 
are exactly constant. Unlike the ordinary GR case, this result 
does not apply to the matter-frame potential O and the den- 
sity contrast which are not strictly frozen due to the generally 
time-varying perturbations of the scalar field. Given the typ- 
ical magnitudes of scalar field perturbations from numerical 
analysis, however, the impact on the evolution of <t> and 5 will 
be rather small. 

In principle, the analysis presented here can also be per- 
formed for more realistic cosmological models, for instance, 
by including DE effects in terms of a cosmological constant. 
Having said that, however, such studies are likely to involve 
a numerical treatment of background quantities. Similar to 
previous considerations, one might also extend the current pa- 
rameterization to the background level in these cases. More- 
over, note that the basic approach of this work is not limited to 
the framework of TeVeS. Since the occurrence of trivial iden- 
tities of the above kind appears as a generic feature of modi- 
fied gravity theories with extra DOFs, these parameterizations 



TABLE II. Summary of relevant perturbation quantities in the con- 
formal Newtonian gauge. 



Symbol 



Description 



Metric potentials in the matter frame 

Metric potentials in the Einstein frame 

Einstein-frame metric perturbation, got = -Irdit 

Perturbation of the scalar field tp 

Auxiliary perturbation related to the gradient of <p 

Perturbation of the vector field A„ 



"P, <t> 

I 

f 

y 

a Longitudinal component of the vector field perturbation^ 

E Auxiliary perturbation related to F /lv , E is gauge-invariant 

Sp Density perturbation of the matter fluid 

SP Pressure perturbation of the matter fluid 

<5 Density contrast of the matter fluid, 6 = Sp/p 

6 Velocity potential of the matter fluid 

£ Shear of the matter fluid 

C„ Sound speed of the matter fluid, C 2 S = SP/dp 

k Magnitude of the conformal wave vector k, i.e. k = |k| 

B 9 New variable absorbing the scalar field DOF, B^ = T7</? 

B y New variable absorbing the vector field DOF, see Eq. (146) 



1 Note that a is related to £ through Eq. 141) . 



should generally prove useful when exploring the cosmologi- 
cal properties of other proposed modifications including gen- 
eralizations of TeVeS and a variety of tensor- vector models. 

As such parameterizations capture the full nature of mod- 
ifications by construction, they further provide an appropri- 
ate framework to model quantities like the metric potentials, 
the gravitational slip or the growth function in a semi-analytic 
fashion. These models might then be used to constrain mod- 
ified effects with the help of available cosmological data, but 
also form a basis for studying other observable imprints like, 
for instance, the integrated Sachs-Wolfe effect. Supplemen- 
tary to numerical analysis, parameterizations based on trivial 
identities thus provide an interesting tool to approach modi- 
fied gravity theories with extra DOFs. 
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Appendix A: Scalar field evolution during tracking 

In the following, I will assume the generalized potential de- 
fined in Eq. < fT~3T > and adopt the notation and definitions used in 
Ref. 12711 . There it has been found that the scalar field evolves 
during tracking as 



= n + 1 



log a, 



(Al) 
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where 



'i = 



d(f> 

dloga 



(A2) 



is approximately constant. Indeed, following the derivation 
presented in Ref. 12711 . one can show that 



P ( 1 + 3w 



1+01 2yU() ^\ 1 - W j 



(A3) 



where/? = ±1 denotes the sign of the scalar field's time deriva- 
tive, i.e. 



P = sgn0. 



(A4) 



To see that the sign in Eq. ( IA3b is chosen appropriately, one 
uses Eq. ( IA1I ) and finds that 

P = sgn Ui - j = sgn 4>i = sgn/3, (A5) 

where I have assumed that \<pi\ <K 1 for the last equality. Note 
that this is justified because of the requirement po » 1 for 
viable cosmological models. 

Here the right-hand side of Eq. (IA3) deserves special atten- 
tion: When evaluating the square root, one needs to take into 
account that the argument's sign does depend on the actual 
choice of w. Therefore, one has 



1 + 3w 



1 



|l+3w| 
II -wl 



(A6) 



which eventually gives the result in Eq. ( l24t . During tracking, 
the field Jl (see Sec. IIIIBI ) evolves as Jl - 2po(l + e), where 



20!+3(l+w) 
log e oc log a, 



(A7) 



and thus 2<f>\ + 3(1 + w) > emerges as a condition for stable 
tracking. For a universe dominated by a cosmological con- 
stant A, one has w = - 1 and therefore p - 1 . Since the time 
derivative of (f> changes its sign when passing from the matter 
to the A era (resulting in p* momentarily going to zero) 1I29I1 . 
it follows that/? = -1 during matter domination. This result is 
in accordance with previous work [fill 13411 where it has been 
shown that <p decreases with time during the matter era. 

Appendix B: Field equations in conformal Newtonian gauge 

Here I will summarize the resulting TeVeS perturbation 
equations for scalar modes in the conformal Newtonian gauge. 
Furthermore, I shall assume a spatially flat spacetime geome- 
try and introduce the fluid's sound speed C s which is defined 
as the ratio between the fluid's pressure perturbation 5P and 
the corresponding density perturbation 8p, i.e. C 2 - 6P/6p. 
As usual, the equations are expressed in Fourier space using 
the conformal wave vector k in accordance with the coordi- 
nate system specified in Eq. (l3~n >. 



1. Matter fluid equations 

The density contrast for scalar modes in the conformal 
Newtonian gauge evolves as 



8= -(l+w)(k 2 8 + 3&)-3-(c 2 s -w)6, (Bl) 



where the velocity potential 9 obeys 



= --(l -2,w)e+-— —5-— — 6»-/t 2 |z + ^, (B2) 
a 1 + w 1 + w 3 

the quantity S denotes the shear of the matter fluid and k — |k|. 
Note that the equations for perturbations of the matter fluid 
remain unaltered compared to the standard case of GR. 



2. Scalar field equation 



The perturbed scalar field equation yields 

j = - 3-7 + -e~% 2 (tp + da) - 2-e*0(3<l> + k 2 () 
b a \ ) a v ' 

+ SnGae-^p [(l + 3C 2 ) 6+(l+3w){V- 2<p)\ 



and 



<p = -—ae-*y + <p¥. 



(B3) 



(B4) 



Here y denotes the perturbation of an auxiliary field intro- 
duced to split the scalar field equation into two first-order 
equations 13111 . 



3. Vector field equation 

The two first-order equations coming from the perturbed 
vector equation are 



K B \e + -e\ = 8jiGa 2 p(l + w) (l - e~ 40 ) (0 - a) 



and 



a = E + y + U - -) a, 



(B5) 



(B6) 



where the auxiliary scalar mode E is gauge-invariant and re- 
lated to F^y, the field strength tensor of A^ Olll . 



4. Generalized Einstein equations 

From the scalar modes of the perturbed generalized Ein- 
stein equations, one obtains the Hamiltonian constraint 



2A: Z <D - 2e^- 1 30 + k 2 l + 3-T 
b \ b 



+ cur*#y - K B k 2 E = %7iGa 2 p(5 - 2<p) 



(B7) 
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and the momentum constraint equation 



Appendix C: Parameterized perturbation equations 



$ + -T - Jlcfxp = AnG^e-^pil + w)8. 



Finally, the two propagation equations read 



(B8) 



In what follows, I will assume a universe which contains 
pressureless matter only, in accordance with the modified EdS 
cosmology introduced in Sec. |IH CI In this case, the fluid's 
pressure components may be neglected, i.e. w — C s - 2 = 0, 
and the background density evolves as p oc a~ 3 . 



1. General case 



60 + 2k 2 (% - e- 4 ^) + 2e- 4 V$ 
+2- (6& + 3^ + 2k 2 l) + 4^(30 + k 2 l) 

-24nGa 2 e-^p(c 2 s 5 - 2wip) = 



The goal is now to express the metric potentials in terms 
of matter fluid variables, using the closure relations from Sec. 
HVBI together with the perturbation equations. Since there re- 
(B9) main three gravitational fields after applying Eqs. (l44l and 
(|46| |. i.e. the two metric potentials and the longitudinal com- 
ponent a of the vector field perturbation, this requires finding 
three linearly independent equations. For the first equation, I 
eliminate the (conformal) time derivative of <f between Eqs. 
(lB6b and ( IB 10b . which leads to 



and 



O 



(l - e^)E - e^¥ - 40a + e** (- +5<f>\% = 0. (CI) 



<D-¥ + , 



l \l + nl 



87rGa 2 p(l+w)I. (BIO) 



Differentiating the above and substituting all remaining time 
derivatives by a suitable combination of the perturbation equa- 
tions eventually gives 



l ^k 2 l + 
B v h 



aT\2 



(1-T(* 



K B 



2 

2 



1 -e 4 * 



+ 87rGa V 40 p 



1 +2- 



l+e 4 M 



l_ e 40 



^2 



-40 



1 + e 4 M - 

a + 40|-l + = lO 

1 - e 4 



„ 4 <* 



4 =(/>[- +5<f>\+- -2-.- 4<f>-+5(f> +5d> 

1 _ ^40 \a la a 1 a 



l-3e 



l_ e 40 \ a 

-B w \lJI(p b 



£ + A7iGa z e-^p 



1 -^T, 



-B (l-e 4 <T 



5 V 



40 "» 

By lyB^ 



1-^_^[ 4 1 



#fl Biz 



1 -e 4< A 



%i + (4 + ^] + (l-^J 



& 



¥ = 0. 



(C2) 



Remembering that £ is related to a through Eq. fiTl i. Eq. ( IC2b obviously connects the perturbations quantities r ¥, <t> and a to the 
matter perturbation variables (in this case the velocity potential only) and the functions B^ and B y . To find a second equation, 
one may start from Eq. ( |46l . Similar as before, I take its time derivative and use the perturbation equations to recast the resulting 
expression into a more convenient form. A bit of algebra then reveals 



B, 



2U B 



-45 H 



-^W + ^l + -^^rt + e - 4 ^B r ^ + 3^-^-^-3^- + - 



2U 



By 



V 



B„b U 



By 



By U 



Byb U 



x [k 2 l + \2nGa 2 e-%Q) + 3 (- - (4 + B 9 )%- ^]- + —Ga 2 e-%(B y - 2?A - 6- 



\Bj b}\ By^ By U 7 Byb U 



u 

~U 



¥ + ^B y Ga 2 e~%5 = 0. 



(C3) 



Finally, the last equation is obtained from eliminating E between Eqs. dCU and (IB7I ). Together with the relations presented in 
Sec. IIVBI one eventually ends up with 



Kb J\k 2 ®- K B ^^k 2 V + K B I "- - 4 

A<t> b 



- %nGa 2 p5 - 2e 4 



1 

BtpBy \ By 



U\b 



By) b" 



1 _ e 40 

8;r 



k> a _ 2e ^ I \ + — t) {k 2 l + \2nGa 2 e-%Q) 



B, 



Ga'e-^p 



<p 



(C4) 



¥ = 0. 
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As is obvious, the set of Eqs. ( IC2b . ( IC3b and ( IC4I > forms a 
closed linear system for the fields ty, c£> and a. Thus the cor- 
responding solution will give the fields as expressions of the 
parameter functions and matter fluid variables only. 



2. Superhorizon limit 

For scales well outside the horizon, the perturbation equa- 
tions can be considerably simplified by considering the limit 
k — > 0, which allows one to begin with Eq. ( f60b . Taking its 
time derivative, the first step is to eliminate y and 6 from the 
resulting equation. This is achieved with the help of 



-(=-*♦£-£ ^^-"M- < C5 » 



u 



6U 



which follows from Eq. (1461 1. and the relation Eq. (1611 . re- 
spectively. Then, making use of Eqs. d44l) and ( f46b together 
with an appropriate combination of the perturbation equa- 
tions, one can derive an equation governing the evolution of 
the Einstein-frame potential <t>. Defining the auxiliary quanti- 
ties 



Ps3---47rGaV 4 *p 
b z 



(C6) 



and 



Q = u 



b b 

the result can be expressed as 



b b 2 , \b- B^b 



B^b 



(C7) 



(B r P + 3Q)® 



+ 



rlf + *H 



b _-' 
'b 



^iBylb -' 



-3 G ^-^ + ^ 



> - B r 
b v B y 



(C8) 



4> = 0. 



As is obvious, an immediate solution to the above is given by 
O = const. Following a similar procedure, it is also possible 
to derive an equation for the scalar field perturbation <p, 



B,, 



_- b\B y ^ lb - U BA- » P 
fl<f>+-\ -f _ 2- + + — + — - 



b U 



b / v- f> /_- M B r 

2 - h + ( l - B *)4> + n-W + bU 



b 
u_ 
u 



U B u 



u 



if + if — 0, 



(C9) 



which takes the general form of a damped harmonic oscillator 
and is further discussed in Sec. I1VDI 
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